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ABSTRACT 

Context. The activity indicators of M dwarfs are distinctly different for early and late types. The coronae of early M dwarfs display 
high X-ray luminosities and temperatures, a pronounced inverse PIP effect, and frequent flaring to the extent that no quiescent level 
can be defined in many cases. For late M dwarfs, fewer but more violent flares have been observed, and the quiescent X-ray luminosity 
is much lower. 

Aims. To probe the relationship between coronal properties with spectral type of active M dwarfs, we analyze the M3.5 and M4.5 
components of the EQ Peg binary system in comparison with other active M dwarfs of spectral types MO. 5 to M5.5. 
Methods. We investigate the timing behavior of both components of the EQ Peg system, reconstruct their differential emission 
measure, and investigate the coronal abundance ratios based on emission-measure independent line ratios from their Chandra HETGS 
spectra. Finally we test for density variations in different states of activity. 

Results. The X-ray luminosity of EQ Peg A (M3.5) is by a factor of 6-10 brighter than that of EQ Peg B (M4.5). Like most other 
active M dwarfs, the EQ Peg system shows an inverse FIP effect. The abundances of both components are consistent within the 
errors; however, there seems to be a tendency toward the inverse FIP effect being less pronounced in the less active EQ Peg B when 
comparing the quiescent state of the two stars. This trend is supported by our comparison with other M dwarfs. 
Conclusions. As the X-ray luminosity decreases with later spectral type, so do coronal temperatures and flare rate. The amplitude of 
the observed abundance anomalies, i. e. the inverse FIP effect, declines; however, clear deviations from solar abundances remain. 
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1. Introduction 

M dwarfs populate the low-mass and low-temperature end of the 
main sequence, with masses of approximately 0.6 M0 and an 
effective temperature of about 3800 K at a spectral type of MO, 

] down to 0.1 Mq and 2200 K at spectral type M9. Conspicuous 
changes in the optical spectra become apparent in the later sub- 
classes as more and more molecules form with decreasing pho- 
tospheric temperature, and molecular bands start to dominate. 

• Dust formation marks the transition from late M dwarfs to the 

' brown dwarf regime of L and T dwarfs. 

Somewhat surprisingly, all M dwarfs seem to show at least 
some activity phenomena, and in particular active M dwarfs of 
spectral types M0-M4 are strong coronal X-ray sources with X- 
ray lu minosities often close to a saturation limit o f log Lx/L^oi ~ 
-3.3 dFleming et alj|l993l; [Pizzolato et al.ll2003l) . Frequent and 
strong flaring with increases in X-ray luminosity by more than 
two orders of magnitude ( e. g. iFavata et al] 120001 : iGtidel et al] 
l2002al : ISchmitt et alj l2008l) is widespread among these stars, 
with additional variability on all timescales and amplitude lev- 
els. Their e mission measur e distributions peak around 6-8 MK 
jRobrade & Schmittll2005l) . thus their average coronal temper- 
atures ar e much higher than solar (peak temperature 1-2 MK, 
see e. g. iBrosius et alj Il996h . Temperatures in excess of sev- 
eral 10 MK are typically measured during larger flare events 
(Robrade & Schmitt 2003). 
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On the other hand, late M dwarfs with spectral types M6- 
M9 are difficult to detect in X-rays at least during quiescence, 
yet their ability to produce transient X-ray luminosity enhance- 
ments by orders of magnitud e during flare s (iRutledge et^aP 
SoOOt ISchmitt & Liefeg I200I IStelzer et al.1 l2006h remains. 
Intermediate objects with a spectral type of ^ M5 have lower 
coronal temperatures than the more active early M dwarfs, the 
permanent state of variability is replaced by periods of quies- 
cence inteiTupted by individual flares, and the overall X-ray lu- 
minosity starts to decrease by more than an order of magnitude at 
aroun d spectral type M5 dGtidel et alj|2002al; iFuhrmeister et al.l 
120071) . 

The abundance patterns observed in stellar coronae also ap- 
pear to be related to changes in the coronal activity level. On the 
one hand, there is the so-called FIP effect observed in (co mpara- 
bly) in active stars like the S un dFeldman & Lamind2000l) or e. g. 
a Cen dRaassen et al.ll200"3l) . with the coronal abundances of ele- 
ments with low first ionization potential (FIP) enhanced with re- 
spect to solar photospheric values. On the other hand, very active 
stars tend to show the inverse effect, i. e. elements with low FIP 
depleted and elements with a high FIP enhanced. i Giidel et"al] 
(2OO2H) andfrefleschi et al. (2005) found that the solar-like FIP 
effect to first disappear and then reverse with increasing activ- 
ity for solar analogs at different ages. This coiTelation can be 
expected to also apply to the sequence of M dwarfs: While 
the active early M dwarfs are well-known to show the inverse 
FIP ef fect (see e. g. Robrade & Schmitt 2005), Fuhrmeiste r et al.l 
d2007l) found the M5.5 star CN Leo to show much less pro- 
nounced abundance anomalies. Unfortunately, high-resolution 
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X-ray spectra for intermediate objects are rare; and due to their 
low X-ray luminosity, it is in most cases impossible to obtain 
spectra suitable for a reasonable abundance analysis for late M 
dwarfs to confirm this trend. 

Stellar coronal X-ray emission is associated with the produc- 
tion of magnetic field s driven by a dynamo mechanism ( Parker 
119551: iMoffattI Il978h . which is thought to be located in the 
transition layers between the rigidly rotating radiative core and 
the outer convection zone that rotates differentially for solar- 
like stars. However, models of the stellar interior predict that 
M d warfs of spectral types la ter than x M3 become fully convec- 
tive ( Chabrie r & Baraffall99 7). leaving no shear layers between 
a radiative core and the convection zone behind. The solar-like 
aw-dynamo thus does not work for these stars, and no dynamo- 
induced X-ray emission would be expected. Yet, the observed 
decline in X-r ay luminosity around s pectra l type M5-M6 fits 
the finding of M ullan & MacDonaldl (1200 1) that the influence 
of magnetic fields can shift the fully convective cut-off towards 
later spectral types. Fully convective stars produce less quies- 
cent X-ray emission, but other activity indicators like Ha emis- 
sion even increase in strength towards the low-mass L dwarfs. 
Alternative dynamo mechanisms, driving the requ ired magnetic 
fields, have been proposed for these objects (e. g. i Durnev et al.l 
Il993h . In these scenarios with or turbulent dynamos, the mag- 
netic fields are sustained by small-scale velocity fields, while 
large-scale rotation plays only a minor role. As a consequence, 
no activity cycles and magnetic structures of small size and with 
a more uniform surface distribution are expected. However, most 
recent non-linear 3 simulations of dynamo action in fully con- 
vective stars (e. g. IChabrier & Kukeiii2006l; iDobler et al.ll2006l; 
[Browning 2008) have shown that also large-scale fields could be 
produced in fast-rotating late M dwarfs, but these models dis- 
agree about whether or not the magnetic fields have axisymmet- 
ric components and/or a cyclic behavior. 

With the required change of the dominating dynamo mecha- 
nism operating in early and late M dwarfs in mind, changes in the 
abundance patterns between these objects are of particular inter- 
est to characterize the coronae induced by the different dynamos. 
The above considerations concerning luminosity demand inves- 
tigations of coronal properties of M dwarfs located close to the 
expected transition. The binary system EQ Peg consists of an 
M3.5 and M4.5 star and thus provides the rare opportunity to 
compare two such stars, which should otherwise be very similar 
since they are co-eval. 

In the following, we present our Chandra HETGS observa- 
tions of EQ Peg, that allows us to perform a separate spectral 
analysis of the two components. We also compare our observa- 
tions of EQ Peg A and B to HETGS spectra of other M dwarfs 
covering a range of spectral types. We propose a classification of 
the coronal properties of active M dwarfs based on their spectral 
type. 

2. The EQ Peg binary system 

The visual binary system EQ Peg consists of two M dwarfs 
at a distance of 6.25 pc with an angular separation of » 6", 
i.e. EQ Peg A with spectral type M3.5 and an apparent mag- 
nitude Va - 10.35 and EQ Peg B with spectral type M4.5 and 
magnitude Vb - 12.4. Both components are we ll -known (opti- 
cal) fla re stars dPettersenll 19761; iLacv et al.lll976h . iJackson et alj 
(11989^ observed microwave emission during quiescence, which 
they attributed to the brighter A component. Norton et al. ( 2007) 
find a photometric period of 1 .0664 days for the EQ Peg sys- 
tem from data of the SuperWASP transit survey, however, due 



Table 1. Observation log for EQ Peg. 



ObsID 


observation time 




exp. time [ ks ] 


8453 


2006- 1 1-28T0 1 :5 1-2006- 1 1-28T22 


48 


29.95 


8484 


2006- 11 -29T 1 6: 33-2006- 11 -29T08 


59 


10.95 


8485 


2006- 11 -30T 1 6 :59--2006- 1 2-0 1 T08 


54 


30.94 


8486 


2006- 1 2-03X07:55-2006- 1 2-03T09 


02 


27.95 



to the large aperture size of SuperWASP, this period cannot be 
clearly assigned to one of the two components of EQ Peg, nor 
to two other stars in the field of view. However, Zboril & Byrr3 
find a r o tation al velocity of 14 km/s for EQ Peg A, while 
Delfosse et alj ( Il998l) analyze only the B component and find 
24.2+ 1 .4 km/s; both v sin / values are rather high and thus con- 
sistent with a short rotation period. 

The EQ Peg system is a strong X-ray and EUV source with a 
number of flares recorded, but pr evious missions were no t able to 
angularly resolve the two stars. Pallavicini et al.l (1 19901) discuss 
two flares observed with EXOSAT. The first one, with an atypi- 
cally shaped lightcurve, was observed bv lHaisch et al.l ( [19^871) in 
the context of a simultaneous EXOSAT and l UE campaign. A 
second large amplitude flare was observed by iPallavicini et al.l 
(Il986h . iKatsova et all (|2002|) observed another large flare with 
ROSAT with simultaneous optical photo metry. All of these flares 
were attributed to the A component. iMonsignori Fossi et al.l 
(Il995h observed another set of flares on EQ Peg with EUVE and 
derived a differential emission measure distribution, based on 
line flux measurements, with a peak at about 10 MK, which is, 
however, not well constrained. A first approach to separate the 
two components in X -rays was undertaken with XMM-Newton 
dRobrade et alj |2004|) . Although the two stars show consider- 
able overlap owing to the instrumental point spread function, 
it was possible to attribute about three quarters of the overall 
X-ray flux to the A component. A subsequent detailed spec- 
tral anal^^sis_W2thou^es^^ the binary has been performed 
bv lRobrade&Schmittl (12005'). Our Chandra HETGS observa- 
tions are the first X-ray observations that allow an unambiguous 
spectral separation of the two binary components. 

3. Chandra HETGS observations of EQ Peg A and B 

3.1. Observations and data analysis 

EQ Peg was observed on four occasions between 2006 
November 28* and December 3"* with a total exposure time 
of 100 ks with the Chandra High Energy Transmission Grating 
Spectrometer (HETGS); a journal of the observations is pre- 
sented in Table [1] The HETGS consists of two grating arms, the 
Medium Energy Grating (MEG) and the High Energy Grating 
(HEG), which provide X-ray spectra in the wavelength range 
between * 1.5 and 30 A with medium (MEG: 0.023 A FWHM) 
and high (HEG: 0.012 A FWHM) resolutionQFor data reduction 
we used the CIAO software version 3.4. 

The separation of the 0* order images of EQ Peg A and B is 
5.8" in the ACIS-S image, and therefore the dispersed spectra of 
the two stars can be easily separated with only a small overlap 
of the HEG and MEG grating arms at the shortest wavelengths 
(< 1.5 A). We extracted P' order HEG and MEG spectra for 



' For a detailed description of the Chandra HETGS 
see the Chandra Proposers' Observatory Guide at 
http://asc.harvard.edu/proposer/POG/html/| 
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each component from each of the four datasets using the stan- 
dard CIAO tools to obtain grating spectra for muhiple sources. 

For our final analysis we created separate HEG and MEG 
spectra for EQ Peg A and B, where we combined all four obser- 
vations and the plus- and minus dispersion orders. Only for EQ 
Peg A have we also inspected the spectra from the individual 
observations. 

For our final analysis, we created separate HEG and MEG 
spectra for EQ Peg A and B, where we combined all four ob- 
servations and the plus and minus dispersion orders. Only for 
EQ Peg A we have also inspected the spectra from the individ- 
ual observations. All line fluxes were measured with our CORA 
program (Ness & Wichmann.2002) . fitting the fines with Moffat 
line profiles given by 

/(^)=w(i+(^)y' (1) 

i. e., modified Lorentzians, with an exponent j3 - 2.4 and a fixed 
line width of 0.02 A for MEG and 0.01 A for HEG spectra. 

3.2. Lightcurves 

We first extracted background- subtracted lightcurves for EQ Peg 
A and B from the dispersed and non-dispersed (0* order) pho- 
tons. In order to obtain equivalent lightcurves of the two stars 
with maximized signal-to-noise, we use total (i. e. 0* order -H 
dispersed) photons of EQ Peg B and the dispersed photons of 
EQ Peg A, scaled by the ratio of total and dispersed photons 
from EQ Peg B (assuming a similar spectral energy distribution 
for the two stars), as the the 0* order of EQ Peg A is strongly 
affected by pileup. The resulting lightcurves are shown in Fig.[T] 
for the four observation intervals. 

EQ Peg A shows several flares. The most conspicuous event 
started at 0:00 UT on December P' with a relatively slow rise 
and lasted about 2 hours. Unfortunately, the flare decay phase 
was not observed completely, since the observation ended be- 
fore quiescence was reached again. At flare peak, the flux in- 
creases by a factor of six, and the flare is strong enough to al- 
low the extraction of separate MEG spectra. Additional smaller 
flares occurred on November 28* at 5:00 UT, with a flux in- 
crease of about a factor of four, lasting about half an hour, and 
on December 3"* at about 9:00 UT, with an increase by a factor 
of 2, and a longer decay of more than 2 hours. The underlying 
quiescent level is difficult to determine owing to a large num- 
ber of small-scale variations and apparently different quiescent 
emission levels during the four individual observations, with the 
base level for the A component at ~ 0.3-0.4 cts/s. Apart from 
pure statistical noise, such low-level amplitude variations can 
be caused either by small-scale flares or the evolution of active 
regions on the stellar surface. On November 29* the emission 
level was much higher than this quiescent level despite a lack of 
clear flare-like signatures, which might be explained by active 
regions newly emerged or rotated into the field of view, while 
on November 30* a slow decrease in count rate is observed that 
might be attributed to a long-duration flare. The B component at 
~ 0.05 cts/s is much fainter and does not show conspicuous short 
duration flares, except on Dec 3'^''. On November 30* the count 
rate of EQ Peg B slowly decreased during the whole 8 hours of 
the observation and again we may have witnessed the decay of a 
flare. 

We computed median and mean count rates in order to com- 
pare the flux levels of the two components outside obvious flar- 
ing periods, and find the quiescence level of EQ Peg B lower than 



that of A by factors of about 6 to 10. Both stars show a decay in 
count rate on November 30*, however, a direct connection of 
flare activity on the two components is difficult to envisage con- 
sidering the distance of ~ 36 AU between A and B, although 
a similar correlated behavior of the system wit h simultaneous 
quiescent and ffaring states has been observed bv lRobrade et al.l 
(2004) with XMM-Newton. These authors also found a lower 
ffux ratio of 4-5 during quiescence and 2.5-3 during flaring 
states. However, due to the limited spatial resolution of XMM's 
EPIC detectors, the exact brightness ratios are less certain. 

We also extracted the lightcurves in different energy bands to 
obtain time-resolved hardness ratios. These give us the tempera- 
ture evolution of the plasma. While there is no unique tempera- 
ture that describes the spectrum (see Sect. I4.2l i. the temperature 
distribution cannot be mapped as a function of time, and the dif- 
ferential emission measure distribution can only by determined 
as an average over the entire observation. Because of pileup in 
the 0* order of the A component, we extracted the hardness 
ratios only from the dispersed photons for both stars. With the 
hard and soft bands chosen to range from 1 .0-10.0 keV and 0.3- 
1.0 keV respectively, the two components show consistent hard- 
ness ratios during quiescence, with values ranging from 0.25 to 
0.5, indicating a similar spectral energy distribution. Averaged 
over the total observing time, the hardness ratio for both stars is 
0.43; a spectral hardening to values up to 0.8 is apparent during 
the major flares on the A component. For our spectral analysis, 
we considered only the total spectrum of EQ Peg B, while for 
EQ Peg A we considered the total spectrum and a "quiescent" 
spectrum, where we excluded the two larger flares on November 
28* and December V\ 

We computed overall X-ray luminosities log Lx of 28.71 for 
EQ Peg A (28.67 for the quiescent state only) and 27.89 for 
EQ Peg B from the integrated dispersed photons in the MEG 
spectra over a wavelength range from 2 to 25 A, i. e. in this spec- 
tral range, the A component is by a factor of 6-7 more luminous 
than th e B component. However, following lKenvon & HartmannI 
(Il995h . the bolometric luminosities of the two stars are log Lboi ~ 
31.90 and 31.28. This results in logLx/i'boi values of -3.19 
(-3.23 during quiescence) and -3.39 for EQ Peg A and B re- 
spectively, i. e. both stars are located in the saturation regime of 
X-ray emission (iPizzolato et al.ll2003l) . 

3.3. Grating spectra 

In Fig. |2] we show the Chandra MEG spectra of EQ Peg A 
(top) and B (bottom). Both stars show well-pronounced con- 
tinua between 5 A and 15 A, indicating larger amounts of high- 
temperature plasma in their coronae. The strongest lines in the 
total MEG count spectra are Ne x and O viii Ly a, other promi- 
nent lines originate from Si xiv and xin, Ne ix, and Fe xvii. The 
iron lines tend to be somewhat stronger in EQ Peg B. The O vii 
triplet lines are clearly detected but suffer from the low effec- 
tive area of the MEG at long wavelengths. The weak O vii ls-3p 
line at 18.63 A is also detected in both stars. After conversion of 
counts to photon fluxes, the O viii Ly a line is the strongest line 
in the spectra, and also the O vii flux is quite high, indicating the 
presence of considerable amounts of low-temperature plasma. 

In order to probe the shape of the continua and to directly 
compare the two binary components, we computed the cumula- 
tive count spectra from 4 A to 26 A from the combined count 
spectra of EQ Peg A and B (top panel of Figure [3j. A steep 
rise of the cumulative spectrum indicates a harder spectrum and 
thus higher temperatures, as observed for EQ Peg A compared 
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EQ Peg ACIS-S pile-up corrected lightcurves 
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Background-subtracted and pileup-corrected lightcurves for EQ Peg A (top) and B (bottom). Binsize for both sources is 



Combined CInandra MEG spectra of EQ Peg A and B 
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Fig. 2. Chandra MEG spectra of EQ Peg A (top, red) and B (bot- 
tom, blue). Both spectra have been rebinned. 



to EQ Peg B. Since this may reflect the strong flares observed 
on the A component, we considered also the quiescent spectrum 
for EQ Peg A. As can be seen from Figure[3](top), this spectrum 
is indeed closer to the B component, however, EQ Peg A still 
appears somewhat harder. When comparing the total spectra of 
A and B, these differences are relatively small, and a KS test re- 
turns probabilities of around 3% that the two stars are different 
in their spectral shape. We note that the presence of strong emis- 
sion lines, e. g., Nex at 12.12 A, leads to steps in the cumulative 
spectrum. However, the depths of these steps seem similar for 
EQ Peg A and B, and the major differences in the cumulative 
distribution can be attributed to the continuum. 

In Tables|2]and[3]we provide lists of the measured line counts 
and resulting photon fluxes for the combined spectra of EQPeg A 
and B, and for the quiescent state of EQ Peg A, which we used in 
our subsequent analysis; note that a sophisticated procedure has 
been applied to deconvolve the blending of the Ne ix triplet and 
its contamination w ith iron lines, mainly from Fe xvii and Fe xix 
jNess et al.ll2003bi) . In addition to the line fluxes of the three 
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Fig. 3. Cumulative distribution (top) of counts in the MEG spec- 
tra of EQ Peg A (dotted line: all data, solid line: with the large 
flare on December 1" excluded) and B (dashed line); and ratio 
B / A (middle panel: all data, bottom: with the two larger flares 
on EQ Peg A excluded). 



triplet constituents, we fitted global scalin g factors to the line 
fluxes of the contamining ions observed by iNess et all (l2003bl) 
in the well-exposed HETGS spectrum of Capella where the con- 
tamination is very strong. Additionally, the Nex line at 12.13 A 
has been corrected to account for a blend by Fe xvii at 12. 12 A. 
According to CHIANTI 5.2 (Landi et al. 2006), using the ion- 
ization balance from iMazzotta et al.. (,1998.) . the line emissivity 
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Fig. 4. Spectral region around the Nex line at 12.13 A and the 
Fe xvii and Fe xxi lines at 12.26 and 12.28 A in the total spectrum 
of EQPeg A. 

of the contaminating line is comparable to the line at 12.26 A, 
and the ratio of these two lines shows no dependence on elec- 
tron density and only a mild dependence on electron tempera- 
ture. Additionally, both transitions go to the ground state and 
share the same electron configuration in their excited states. We 
thus subtracted the line counts measured in the line at 12.26 A 
from the Ne x line as a proxy for the unresolved blending line 
at 12.12 A. Fig.|4]illustrates the fit for the total HEG and MEG 
spectra of EQ Peg A. 

4. Results 

4. 1 . Coronal densities 

Using the forbidden and intercombination lines of the helium- 
like triplets of silicon, magnesium, neon, and oxygen, we com- 
puted the density-sensitive f/i-ratios of these ions for EQ Peg A 
and B. In Table |4] we list these ratios with the respective peak 
formation temperatures of the ion and the derived electron den- 
sities He for EQ Peg A and B. For the conversion of the measured 
/// ratios to densities, we used the relation 



with the low-density Umit Ro and the critic al density N c wher e 
we adopted the values from Table 5 in 'N ess et all (|2002|) . 
Theoretical modeling of the // / ratios still suffers from the in- 
completeness of current atomic databases, as contributions from 
satellite lines and effects from dielectronic recombination have 
to be taken into account. 

The uncertainty in the density derived from O vii is large, as 
a result of the low effective area at long wavelengths. While the 
Ne IX triplet around 13. 5A is well exposed, severe contamination 
complicates the line flux measurements, especially of the Ne ix i 
Une at 13.54 A dNess e t al. 2003b), see also Sect.|33] The den- 
sity value derived for EQ Peg A is higher than the one derived 
from O VII. Since Ne ix traces higher coronal temperatures, this 
may indicate a higher pressure in the hotter plasma regions, how- 
ever, the uncertainties do not include uncertainties from the de- 
blending procedure. The ///-ratio derived for the B component 
is - nominally - in the low-density limit region, however, the er- 
ror is so large that the same density derived for the A component 



Table 3. Measured line counts and photon fluxes in the Chandra 
MEG spectra of the quiescent state of EQ Peg A. 



^[A] 


ion 


counts 


fluxes 


4.73 


S XVI 


8.6+4.1 


0.25+0.12 


5.04 


Sxv 


22.6±5.8 


0.70+0.18 


6.18 


Si XIV 


143.7+13.2 


1.54+0.14 


6.65 


Sixm 


317.7±19.1 


2.87+0.17 


6.69 


Si XIII 


57.0+10.1 


0.52+0.09 


6.74 


Si XIII 


193.6+15.6 


1.60+0.13 


8.42 


Mgxii 


175.4+15.3 


1.51+0.13 


9.17 


Mgxi 


98.7+11.8 


1.54+0.18 


9.23 


Mgxi 


27.5+8.1 


0.44+0.13 


9.31 


Mgxi 


31.9+8.4 


0.54+0.14 


10.24 


Nex 


155.6+14.0 


2.49+0.22 


11.17 


Fe xxiv 


20.6+7.6 


0.44+0.16 


11.55 


Neix 


90.3+11.5 


2.02+0.26 


11.74 


Fe xxiii 


45.9+9.0 


1.12+0.22 


11.78 


Fe xxii 


42.4+9.0 


1.11+0.24 


12.13 


Nex 


813.8+28.5 


23.75+0.83 


13.45 


Neix 


293.4+17.1 


13.67+0.80 


13.55 


Neix 


61.1+7.8 


2.76+0.35 


13.69 


Neix 


169.1 + 13.0 


8.69+0.67 


14.21 


Fe xviii 


82.7+10.6 


5.53+0.71 


15.01 


Fe XVII 


266.3+16.9 


16.45+1.05 


15.26 


Fe XVII 


92.8+10.7 


7.39+0.77 


16.01 


VIII 


184.8+14.0 


13.69+1.05 


16.78 


Fe XVII 


114.3+11.0 


10.80+1.04 


17.05 


Fe XVII 


139.9+12.6 


14.73+1.33 


17.10 


Fe XVII 


137.2+12.4 


14.70+1.33 


18.63 


Ovii 


30.3+5.6 


4.55+0.87 


18.97 


VIII 


414.0+20.5 


72.49+3.60 


21.60 


Ovii 


64.5+8.3 


23.76+3.06 


21.81 


Ovii 


14.0+4.0 


5.91 + 1.69 


22.10 


Ovii 


24.0+5.1 


11.62+2.47 


24.79 


N VII 


15.0+4.3 


6.47+1.85 



Table 4. /// ratios (in photon units) and coronal densities for 
EQ Peg A and B from total Chandra MEG spectra. 







EQPe 


gA 


EQPe^ 




ion 


logr 


f/i 


log He 


f/i 


log 


Ovii 


6.3 


1.69+0.60 


10.62+"-^^ 


2.81+0.87 


< 10.98 


Neix 


6.6 


2.38+0.32 


11.40t^f 


3.53+1.33 


< 11.51 


Mgxi 


6.8 


1.01+0.35 


12 99+0 27 


1.68+1.04 


< 13.28 


Sixm 


6.95 


3.28+0.62 


< 11.17 


5.08+3.64 


< 13.52 



is also consistent with the data. The Mg and Si triplets are much 
stronger affected by satellite lines. Additionally, the Mg xi triplet 
is blended with higher-order Unes f rom the Lyman series of Ne x 
as discussed by iTesta et alj ( 12004 !). The flux in the magnesium 
triplet lines is too low in both components of the EQ Peg bi- 
nary for a deblending procedure to yield any useful refinements. 
Since neon turns out to be far more abundant than magnesium 
(see Sect. |4.2| i, the blending contributions from neon amounts 
to a significant fraction of the line fluxes, and we consider the 
high densities derived from magnesium not reliable. Consistent 
with that interpretation, the Si xiii triplet lines are consistent with 
being in the low-density limit, which is not surprising since the 
density-sensitive range for this ion is higher than the densities 
normally encountered in stellar coronae. However, the Si triplet 
is also problematic because of a steep increase of the effective 
area in both the HEG and MEG that occurs exactly at the wave- 
length of the forbidden line. 
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Table 2. Measured line counts and photon fluxes (in 10 ^ photons cm ^ s ') in the combined Chandra MEG spectra. Flux predic- 
tions from the reconstructed DEMs (see Setc. 14.21 ) are also listed. 



A 

[A] 


ion 


counts 


EQ Peg 
photon fluxes 


A 

method 1 


method 2 


counts 


EQ Peg 
photon fluxes 


B 

method 1 


method 2 


4.73 


S XVI 


22.3±5.6 


0.61+0.15 


0.42 


0.32 


— 


— 


— 


— 


5.04 


Sxv 


28.7+6.3 


0.84+0.18 


0.92 


0.94 


— 


— 


— 


— 


6.18 


Si XIV 


176.6±14.6 


1.78+0.15 


1.92 


1.86 


30.4+5.9 


0.31+0.06 


0.21 


0.13 


6.65 


Sixin 


321.6±19.7 


2.74+0.17 


2.71 


2.75 


40.4+6.8 


0.34+0.06 


0.36 


0.37 


6.69 


Si XIII 


61.7+10.6 


0.53+0.09 


0.72 


0.73 


5.0+3.4 


0.04+0.03 


0.10 


0.10 


6.74 


Sixm 


202.5+16.2 


1.58+0.13 


1.31 


1.33 


25.7+6.0 


0.20+0.05 


0.19 


0.19 


8.42 


Mgxn 


178.9+15.9 


1.45+0.13 


1.58 


1.60 


30.2+6.0 


0.24+0.05 


0.25 


0.23 


9.17 


Mg XI 


107.7+12.4 


1.58+0.18 


1.54 


1.50 


26.0+5.6 


0.38+0.08 


0.32 


0.34 


9.23 


Mgxi 


33.9+8.2 


0.51+0.12 


0.36 


0.35 


6.7+3.4 


0.10+0.05 


0.08 


0.08 


9.31 


Mg XI 


34.6+8.4 


0.55+0.13 


0.81 


0.78 


11.2+3.9 


0.18+0.07 


0.17 


0.18 


10.24 


Nex 


160.8+14.4 


2.42+0.22 


3.05 


3.07 


29.4+5.91 


0.44+0.09 


0.34 


0.36 


11.17 


Fe XXIV 


47.2+9.7 


0.95+0.19 


0.48 


0.44 


— 


— 


0.02 


0.01 


11.55 


Neix 


101.3+12.3 


2.13+0.26 


1.80 


1.84 


19.5+5.12 


0.41+0.11 


0.28 


0.28 


11.74 


Fexxm 


65.1 + 10.5 


1.50+0.24 


0.91 


1.00 


— 


— 


0.13 


0.09 


11.78 


Fexxn 


44.9+9.6 


1.11+0.24 


1.12 


1.29 


— 


— 


0.11 


0.17 


12.13 


Nex 


891.8+29.9 


24.47+0.82 


24.00 


24.06 


104.8+10.6 


2.88+0.29 


2.75 


2.92 


13.45 


Ne IX 


320.2+17.9 


15.89+0.88 


17.02 


16.58 


51.4+7.5 


2.26+0.33 


2.66 


2.56 


13.55 


Neix 


77.2+8.8 


3.28+0.37 


3.67 


3.57 


9.1+3.0 


0.39+0.13 


0.56 


0.56 


13.69 


Neix 


184.0+13.6 


8.90+0.66 


9.53 


9.28 


32.0+5.7 


1.55+0.28 


1.45 


1.46 


14.21 


Fe xvm 


92.2+11.1 


5.80+0.70 


4.77 


4.86 


18.7+4.6 


1.18+0.29 


0.70 


0.76 


15.01 


Fe XVII 


270.1 + 17.3 


15.72+1.01 


18.74 


18.51 


43.2+7.0 


2.51+0.40 


3.29 


3.31 


15.26 


Fe XVII 


102.2+9.3 


6.60+0.60 


5.36 


5.29 


26.5+5.7 


1.71+0.36 


0.95 


0.95 


16.01 


VIII 


193.1 + 14.5 


13.71+1.03 


8.40 


8.33 


27.2+5.7 


1.93+0.40 


0.89 


1.03 


16.78 


Fexvn 


116.4+11.3 


10.36+1.00 


11.41 


11.25 


18.3+4.6 


1.63+0.41 


2.05 


2.04 


17.05 


Fexvn 


144.7+12.7 


14.36+1.26 


14.51 


14.30 


24.8+5.2 


2.46+0.52 


2.60 


2.59 


17.10 


Fexvn 


139.3+12.4 


14.07+1.25 


12.10 


11.93 


26.2+5.4 


2.65+0.55 


2.15 


2.16 


18.63 


Ovn 


33.9+6.2 


4.78+0.87 


3.37 


3.43 


7.3+3.1 


1.03+0.43 


0.51 


0.48 


18.97 


vm 


435.8+21.1 


71.90+3.48 


70.50 


69.90 


63.4+8.2 


10.46+1.36 


7.77 


8.82 


21.60 


Ovn 


69.6+8.6 


24.21+2.98 


30.53 


31.30 


7.6+2.8 


1.56+0.91 


4.78 


4.40 


21.81 


Ovn 


14.2+4.0 


5.64+1.59 


6.83 


7.07 


3.7+2.0 


1.47+0.79 


1.45 


0.99 


22.10 


Ovn 


24.0+5.1 


10.95+2.33 


18.08 


18.70 


10.3+3.5 


4.70+1.6 


3.71 


2.62 


24.79 


N vn 


19.3+4.5 


7.84+1.85 


7.84 


7.84 


2.0+1.7 


0.80+0.72 


0.80 


0.80 



The flux in the Ne ix triplet of EQ Peg A is large enough 
to investigate not only the density in the combined spectrum, but 
also in the individual datasets. In Fig.|5]we show the correspond- 
ing spectral region for each of our four data sets. While differ- 
ences in the strength of the intercombination line at 13.55 A are 
apparent in the individual datasets, no changes in density can 
be established at a statistically significant level. In addition to 
the complete dataset 8485, we show also the spectrum extracted 
for the large flare (shaded in gray in the bottom left panel of 
Fig- El). Dividing this observation in flare and non-flaring states 
gives ///-ratios of 1.13+0.59 and 2.24+0.61, or density values 
log rie of 12.07^Q3j and 11.49;^y43, respectively. For the com- 
bined quiescent spectrum, we find /// - 2.77+0.42, correspond- 
ing to log lie - W. 1 3^o o9> which is lower but still consistent with 
the total spectrum (see Table |4| and clearly lower than during 
the large flare. This indicates activity-related density variations 
as observed by Maggio & Ness (2005) on the active M dwarf 
AD Leo. Unambiguous measurements of coronal density varia- 
tions have so far only been observed during very large fl ares on 
the active M dwarfs Prox ima Centauri and CN Leo (Giid el et alj 
l2002al : ILiefke et aT* 2007). The two density values derived here 
for flare and quiescence are consistent with each other at the 2cr 
level, however, the density values observed during flares on ac- 
tive stars in general - although consistent with the quiescent state 
- tend to result in higher densities (see e. g. .Mitra-Kraev & Ness. 
l2006h . 



EQ Peg A: Ne IX triplet 




13.4 13.5 13.6 13.7 13.8 13.4 13.5 13.6 13.7 13.8 
MA] A [ A ] 



Fig. 5. The Ne ix triplet in the MEG spectra of EQ Peg A during 
the four observations. On observation No. 8485 the flare spec- 
trum only is overlaid in gray. 



We also note the different levels of contamination of the neon 
triplet by Fexvn (13.8 A) and Fexix (13.51 A) in the dififerent 
datasets. During the large flare, this can be attributed on the one 
hand to the higher coronal temperatures increasing the amount 
of Fe XIX compared to other ionization stages of iron, and on the 
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Other hand possibly also to an increased Fe/Ne abundance ratio 
during the flare, as discussed below. 

4.2. Differential emission measure 

In order to derive absolute abundances we investigated the dif- 
ferential emission measure (DEM) of the two stars from the line 
fluxes of the strongest density-insensitive lines listed in Table |2l 
For EQ Peg A, the HEG spectrum can also be used in addition 
to the MEG spectrum, while for EQ Peg B only the MEG spec- 
trum provides data with a reasonable signal-to-noise ratio. In a 
two-step procedure we first reconstructed the DEM utilizing an 
abundance-independen t approach similar to the one proposed by 
ISchmitt & NessI (|2004|) . In a second step we then determined the 
abundances from the ratios of measured line fluxes to those pre- 
dicted by the DEM. 

We used a volume differential emission measure defined as 



DEM(T) = ni 



dV 
df 



(3) 



in units cm"^ K ' . The DEM is related to the theoretical photon 
flux fji (in photons cm"^ s"') of a spectral line originating from 
an atomic transition from the upper level j to the lower level / in 
the ion by 



47rj2 



G(T) DEM(T) dT 



(4) 



where Ax is the abundance of element X relative to solar photo- 
spheric values and the line contribution function G{T) with 



GiT) = Aji 



tijiX^^) n{X^+) n(X) n{H) 1 
n{X^+) niX) niH)^ n, ^ 



(5) 



Aji being the Einstein coefficient for spontaneous emission of 
the transition, n j{X^^) / n{X^^) the level population of the upper 
level, n(X^^)/n(X) the ionization balance, n{X) I niH)^ the solar 
photospheric abundance, and n(H)/ng the proton-to-electron ra- 
tio, which is ^ 0.83 for coronal plasmas where hydrogen and 
helium are usually fully ionized. While the DEM describes the 
properties of the emitting plasma, the line contribution function 
summarizes constants and the u nderlying atomic physics. We 
used C HIANTI 5.2 (jL andi et al. 2006) and the ionization bal- 
ance oflMazzott a et alj ^1998) to calculate G{T) for each line. 

Our abundance-independent reconstruction of the DEM is 
based on line ratios, and we used the ratios of the H-like Lya 
and the He-like resonance lines of oxygen, neon, magnesium, 
and silicon, plus several ratios of adjacent ionization stages of 
iron from Fe xvii to Fe xxiv. Since these line ratios involve only 
lines of the same element, the elemental abundances cancel out, 
and each ratio poses a constraint on the shape of the DEM. 
The choice of these lines yields a temperature coverage from 
^ 2 to 20 MK, however, the low effective area at longer wave- 
lengths with low-temperature lines like O vii and O viii results in 
relatively large uncertainties at temperatures below 5 MK. For 
EQ Peg A, a total of 13 ratios have been used, while some of 
the ratios consisting of higher ionization stages of iron could 
not be used for the B component because of low signal, leav- 
ing only 9 usable ratios. In order to fix the absolute level of 
the DEM in addition to constraining the shape by the ratios, but 
still independently of a reference element like iron, we force the 
DEM and the theoretical continuum emissivity as implemented 
in CHIANTI to reproduce the observed continuum flux. For both 
components of EQ Peg, we used the wavelength region around 



7.5 A, where the effective area of the MEG reaches its maxi- 
mum, and the spectrum is, according to the CHIANTI database, 
essentially line-free, exept for a few weak aluminum and mag- 
nesium lines. We assume that any remaining emission at these 
wavelengths originates from the continuum and neglect the con- 
tribution of any further weak transitions. As a result of merging 
the four datasets, the determined DEMs constitute an average 
over these observations with the different states of activity and 
quiescent and flaring periods of the two stars. 

We used two slightly different approaches to model the 
DEM, that we parameterize by polynomials of different orders. 
In our first method (method 1) we fitted the polynomial param- 
eters as log DEM, and with our second method (method 2) we 
model the linear DEM, both as a function of log T. While for 
method 1 no boundary conditions are needed, for method 2 we 
forced the value of DEM to be zero at two variable temperatures, 
and positive betwee n these temperatures. For a more d etailed de- 
scription we refer to lLiefke & SchmittI din preparation ). 

For EQ Peg A, both methods give acceptable best-fits to 
the selected line ratios assuming 3'^'' and 5* order polynomials, 
yielding values of reduced;t'^ of 2.06 and 2.36 for methods 1 and 
2, respectively. These DEMs are shown in the top panel of Fig.|6l 
The polynomials are an analytical representation of the shape 
of the DEM, which is only valid for certain temperature range 
which is at first restricted by the temperature range covered by 
the emissivities of the lines implemented in the fit. In Fig.|6]we 
additionally plot the polynomials with solid lines in the temper- 
ature range we consider as well-constrained from the available 
lines, and with dashed or dotted lines outside this range. Since 
the high-temperature regime is better represented by the avail- 
able ratios than the low-temperature range (see above), the two 
methods yield best agreement at high temperatures. However, 
the peaks in the DEM disagree from the two methods, and the 
poor constraints from temperatures log T below 6.5 lines make it 
difficult to quantify the exact amount of cool material; this may 
even falsely hint at a DEM peak at a lower temperature. 

For EQ Peg B, this becomes even more apparent, see the 
bottom panel of Fig. |6l The best-fit results are obtained with 4* 
order polynomials with both methods, with reduced values 
of 1.67 and 1.85; however, the two approaches give consistent 
results only in the range of 6.8 < log T < 7.2. The DEM recon- 
structed from method 2 shows an unphysical increase towards 
lower temperatures due to the fact that no ratios are available to 
adequately constrain the low-temperature slope. 

When introducing higher-order polynomials (orders > 5), 
the DEMs of both stars start to develop two-peaked structures 
with maxima at log T » 6.2-6.4 and 6.8-7.2. This would cor- 
respond to a low- and a high temperature component, where the 
first one causes the strong oxygen line fluxes, and the hotter com- 
ponent is related to the typical active M dwarf corona and the ob- 
served flares. However, the DEMs constructed with higher-order 
polynomials are often unstable, and the location and strength of 
certain structures often depends only on a single line flux mea- 
surement. Thus we cannot be certain about the reality of these 
structures, and we prefer the "simpler" DEMs with fewer free 
parameters. 

Next, we used the DEMs from these two approaches for 
the two stars to predict line fluxes. Since the DEMs are deter- 
mined independently of the abundances and are normalized to 
the continuum, any discrepancies between predictions and mea- 
surements of line fluxes must be due to different abundances than 
those assumed in predicti ng the line fluxes from a given DEM 
(here lAsplund et alj|2005h . 
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EQ Peg A: differentiol emission meosure 



E 

" 44.5 




method 1, red. x 2.06 

method 2, red. x' 2.36 



43.51 , 

6.0 



6.4 6.6 6.8 7.0 7.2 

log T [ K ] 

EQ Peg B; differential emission measure 




method 1, red. x 1-67 

method 2, red. 1-85 



6.6 6.8 
log T [ K ] 



Fig. 6. Differential emission measures of EQ Peg A (top) and 
B (bottom) for the total observation time (i. e. flare and quies- 
cence). 



In Table|5]we list our results for the abundances with statisti- 
cal uncertainties derived from the measurement uncertainties of 
the line fluxes; note that the quoted errors only reflect the statisti- 
cal precision of the line flux measurements, but not the accuracy 
of the DEM nor the uncertainties from the continuum flux nor 
any errors arising from the atomic physics parameters used. The 
errors in the continuum amount to 2.5% and 8.2% for EQ Peg A 
and B, respectively, and affect only the absolute abundance level. 
While we experience great uncertainties in determining the ex- 
act shape of the DEM, the determination of abundances seems 
to be only relatively little affected. Polynomials of different or- 
ders always give abundances consistent within the errors and the 
abundances derived with method 1 are always consistent with 
those obtained from method 2. Similarly. iNess & JordanI (l2008l) 
performed a detailed emission measure analysis of e Eridani, 
including ionization stages from many different ions. Close to 
the p eak, their emission measure distribution agrees with that by 
ISanz -Forcada et al. (2004), but is more accurate at low tempera- 
tures. In spite of the differences, both find similar abundances. In 
general, the derived abundances turn out to be quite insensitive 
to changes in the shape of the DEM, polynomials of different 
orders always give abundance value consistent within the errors. 
Additionally, the abundances derived with me thod 1 are al ways 
consistent with those obtained from method 2. lLiefk:e & SchmittI 
din preparation!) assess that the uncertainties on the abundances 
introduced by uncertainties on the shape of the DEM do not ex- 
ceed 5%. 

In Table l5] we show that the absolute abundances relative 



Table 5. Absolute coronal abundances (relative to lAsplund et al.l 
l2005 b of EQ Peg A and B for the total observation time, deter- 
mined with the DEM methods. 





EQPej 


;A 


EQP 


2gB 


element 


method 1 


method 2 


method 1 


method 2 





N 


0.42±0.10 


0.39+0.09 


0.1 6h 


t0.15 


0.20 


+0.18 







0.42+0.04 


0.40+0.04 


0.15H 


t0.09 


0.26 


+0.08 




Ne 


1.14+0.07 


1.07+0.07 


0.93H 


t0.07 


0.83 


+0.07 




Mg 


0.22+0.01 


0.20+0.02 


0.28H 


t0.04 


0.27 


+0.04 




Si 


0.51+0.02 


0.46+0.02 


0.54H 


t0.07 


0.52 


+0.07 




S 


0.64+0.16 


0.58+0.21 












Fe 


0.19+0.02 


0.18+0.01 


0.20H 


t0.02 


0.17 


+0.02 



t dAsplund et all (|2005) of EQ Peg A and B are roughly con- 
sistent for all elements, yielding a low metallicity compared to 



solar This is typic al for active M dwarfs (co mpare e. g. with the 
values derived bv lRobrade & SchmittI '2005, from global fitting 
approaches to CCD X-ray spectra). Nitrogen, oxygen, and neon 
may be somewhat more abundant in the A component, however, 
not at a statistically significant level; also note that the uncer- 
tainties for N and O are quite large. While the high-FIP element 
neon is about solar, the low-FIP element iron shows the lowest 
abundance value, and the other elements also follow the typical 
inverse FIP effect. 

The computed abundances and DEMs can now be used to 
model the spectra of EQ Peg A and B. We list the calculated the- 
oretical line fluxes together with the measured fluxes in Table|2] 

4.3. Emission measure-independent abundances 

Given the ill-posed nature of the DEM determinations and the 
resulting uncertainties in the DEM distributions, we also investi- 
gated emission measure-independent abundance determinations. 
The strongest lines in stellar coronal X-ray spectra (and thus 
easy to measure even in spectra with low signal) originate from 
the H-like and He-like lines of carbon, nitrogen, oxygen, neon, 
magnesium, and silicon, as well as Ne-like Fexvii. In Fig. [T] we 
show the line emissivities for these ions as a function of tem- 
perature. While these lines are formed over a broad tempera- 
ture range, their emissivity peaks are relatively narrow. Those 
lines for which the line contribution functions have similar peak 
formation temperatures yield a similar temperature dependence. 
Ratios of their line fluxes depend only weakly on temperature, 
but stron gly on the a b undan ce ratio of the respective elements. 
Akeady lActon etaTI (1 1975b proposed to determine the solar 
coronal Ne/O abundance ratio from the ratio of the measured 
fluxes of the Ne ix resonance and O viii Ly a lines. This method 
yields a comparably high precision (because of smaller measure- 
ment uncertainties), and also a higher accuracy (because no sys- 
tematic uncertainties from the DEM are introdu ced). Recent ex- 
amples of such ap proaches were presented by ISchmitt & NessI 
(2003) and|Drakel& Testa (2005). 

While for a given element, the peak formation temperature of 
the H-like ion always exceeds that of the corresponding He-like 
ion, the peak formation temperatures of the H-like and He-like 
ions shift to higher temperature with increasing atomic mass. 
Therefore, pairs of H-like and He-like ions of elements with dif- 
ferent atomic mass yield similar peak formation temperatures, 
and their contribution functions have similar shapes. Further ex- 
amples for such pairs are the H-like lines and He-like lines of 
carbon and nitrogen, nitrogen and oxygen, neon and magne- 
sium, and magnesium and silicon respectively (Parkinson 19771). 
Ho wever, these ratios are n ot perfectly temperature-independent, 
and lDrake & Testal (l2005h have refined Acton's method by com- 
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Table 6. Coefficients for Eqns |6] and [T] to obtain temperature- 
independent abundance ratios. 



6.5 7.0 
log T [ K ] 



Fig. 7. Line emissivity curves for different H-like and He-like 
ions and Fexvii at 15.01 A from C HIANTI 5.2 (iLandiet alJ 
120061) with the ionization balance from lMazzotta et alj (Il998l) " 



puting a linear combination of Ne x Ly o- and the Ne ix r line in 
order to construct a new contribution function of Ne lines that is 
more similar in shape to the O viii Ly a line. 

This suggests to ask in general, which linear combinations 
of strong lines without significant blends (i. e. hydrogen- and 
helium-like lines and the strongest lines of Fe xvii as mentioned 
above) of two certain elements yield the smallest temperature 
residuals in their normalized contribution functions. We thus 
write for the abundance ratio of two elements with respective 
abundances A 1 andA2 



A2 



/(Af ) + Ci-/(Af) 

./(A|)' 



C2-/(Af-i) + C3 



(6) 



) denoting the measured line fluxes / 



with /(Af-i) and /(Af- 
of the He-like resonance lines and /(A'^), and /(A|) the corre- 
sponding H-like Lya lines of two elements Ai and A2, respec- 
tively, and 



A_ ^ C4-(/(A^-') + C5-/(A^)) 
Ape 2 ./(Fe XVII) 



(7) 



where ^/(Fexvii) corresponds to /(Fexvii 15.01 A) -H 
/(Fexvii 16.78 A) + /(Fexvii 17.05 A) + /(Fexvii 17.09 A), 
to determine relative abundance ratios. 

To obtain the coefficients Ci, C2, C3, C4, and C5, we per- 
formed a minimization of the temperature residuals of the corre- 
sponding linear combination of the theoretical emissivities e for 
the involved lines over a given temperature range: 

= ^(6,(Af-', r,.)+Cre,(Af , r,.)-C2-e,(Af-', r,.)-C3-e,(Af , T, 

(8) 

for the coefficients in Eqn|6l and 

= ^(e,(A^"', 7^/) + C5 ■ e,(AZ, T,) - i- • ^eKFe xvii, T^f 

(9) 

for the coefficients in Eqn|2l 

In Table|6]we list these coefficients for the temperature range 
where the line emissivities of all involved lines exceed 1% of 



ratio 


Ci 


C2 


C3 


C4 


Cs 


N/C 


+0.13 


-0.07 


+0.73 






0/N 


+0.30 


+0.01 


+0.93 






Ne/O 


+0.02 


-0.17 


+0.69 






Mg/Ne 


+0.18 


-0.08 


+0.87 






Si/Mg 


+0.32 


+0.05 


+0.86 






S/Si 


+0.42 


+0.15 


+0.85 






Ne/Fe 








+34.71 


+0.46 


Mg/Fe 








+67.73 


-0.30 



their peak emissivity. In the given form, the coefficients con- 
vert the measured line fluxes in photon flux units to absolute 
abundance ratios, i. e. independent from any set of solar photo- 
sph eric abundances; w e used line emissivities from CHIANTI 
5.2 (lLandietal.1120061) . 

There are still residuals in temperature for these linear com- 
binations, and their amplitudes differ for each ratio. We found 
the lowest residuals for the Si/Mg ratio and the largest ones for 
the Ne/O ratio. In Fig. [8] we illustrate the best-fit linear combi- 
nation of emissivities that yield the Si/Mg abundance ratio, and 
for the Ne/O ratio we refer to Fig. 3 of Liefke & Schmitt (2006|), 
who had used the same coefficients as listed here to calculate the 
Ne/O abundance ratio of a Centauri. The first six ratios involve 
only H-like Ly a and He-like resonance lines. However, the S/Si 
ratio should only be used when the sulfur lines can be measured 
with reasonable signal-to-noise in Chandra HETGS spectra. 
Fe XVII can be matched with linear combinations of hydrogen- 
and helium-like lines of neon or magnesium. The residuals are 
much smaller for Mg, however, the Mg lines are often weak, in- 
troducing large measurement uncertainties. Especially the weak 
Mg XI line carries a large weight, and for the determination of Fe 
abundances we prefer the Ne/Fe ratio despite its larger residuals. 

The coefficients of the linear combinations as listed in 
Table |6]can be adjusted in some cases where certain line fluxes 
are not available (e. g. the C v lines can only be measured with 
the Chandra LETGS), or are extremely weak (e. g. Mg xi or N vi, 
depending on temperature). The line fluxes /(A^"') and /(A|) 
with Z(Ai) >Z(A2) provide the major constituents of the linear 
combinations for the ratios not involving Fe. 

For example, the N vi resonance line is often weak in stars 
with hot coronae and is difficult to measure in Chandra MEG 
spectra, however, this line contributes much less to the nitrogen 
linear combination in the O/N ratio than the N vii Ly a line (i. e. 
C2 is much smaller than C3). We readjusted C\ , and C3 assuming 
the fixed value C2 = and found only marginal changes for the 
values of the other coefficients listed in Table |6] Although the 
value of Ci is also quite low for the Ne x Ly a line in the Ne/O 
jatio, the situation is somewhat different. If we readjust the other 
)) coefficient after setting Ci = 0, the residuals are much larger, 
especially at high temperatures where the emission measure dis- 
tributions of active stars usually have their maximum. 

Finally, the coefficients can be re-computed for the sum of 
the He-like resonance, intercombination, and forbidden lines for 
cases where these lines can not be resolved (e. g. Si xiii with the 
Chandra LETGS). 

We applied Eqns |6] and [T] with the coefficients from 
Table [6] multiplied by the solar photospheric abundances of 
Asplund et al.l (l2005l) . to the line fluxes measured from the spec- 
tra of EQ Peg A (including and excluding periods of flare activ- 
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Linear combinations of contribution functions: Si and Mg 




log T [ K ] 



Table 8. Abundance ratios relative to lAsplund et alj (l2005h of 
EQ Peg A and B, determined from the absolute abundances 
listed in Table |5] 





EQPe 


gA 




EQ Peg B 




ratio 


method 1 


method 2 


method 1 


method 2 


O/N 


1.00+0.25 


1.03+0.26 


0.94H 


tO.94 


1.30 


+ 1.24 


Ne/O 


2.71+0.31 


2.68+0.32 


6.20H 


t3.75 


3.19 


+0.94 


Mg/Ne 


0.19+0.01 


0.19+0.02 


0.30H 


t0.05 


0.33 


+0.06 


Si/Mg 


2.32+0.14 


2.30+0.25 


1.93H 


tO.37 


1.93 


+0.39 


S/Si 


1.25+0.32 


1.26+0.46 










Ne/Fe 


6.00+0.73 


5.94+0.51 


4.65H 


tO.58 


4.88 


+0.71 


Mg/Fe 


1.15+0.13 


1.11+0.13 


1.40H 


tO.24 


1.59 


+0.30 



Fig. 8. Best-fit linear combination of the contribution functions 
of the H-like Ly a and He-like resonance lines of silicon and 
magnesium. 

Table 7. Abundance ratios (relative to lAsplund et al ."2005') ap- 
plying temperature structure-independent linear combinations of 
coronal emission lines. Values given in brackets for EQ Peg A 
correspond to the quiescent state. 



ratio 


EQ Peg A 


EQ Peg B 


O/N 


0.83+0.20 


(1.00+0.29) 


0.81+0.72 


Ne/O 


2.48+0.31 


(2.02+0.16) 


2.18+0.41 


Mg/Ne 


0.20+0.02 


(0.19+0.02) 


0.32+0.07 


Si/Mg 


2.52+0.21 


(2.56+0.25) 


2.05+0.44 


S/Si 


1.26+0.20 






Ne/Fe 


7.30+ 0.64 


(5.99+ 0.34) 


5.45 + 0.78 


Mg/Fe 


1.27+ 0.17 


(1.09+ 0.19) 


1.63+0.43 



ity) and of EQ Peg B and list the resulting abundance ratios in 
Table 13 

Within the uncertainties, the derived abundance ratios are 
consistent for all three datasets. Only the Ne/Fe ratio is con- 
spicuously lower in the quiescent spectrum of EQ Peg A com- 
pared to the total spectrum, and again lower for the more inactive 
EQ Peg B. Flare- or activity -related changes in the abundance ra- 
tios cannot be determined from count statistics; however, if real, 
the differences may be caused by a decreased iron abundance 
level or a higher neon abundance on the more active A compo- 
nent, especially during flares. 

In Table[8]we list the abundance ratios formed from the abso- 
lute values listed in Table|5]for comparison. They are consistent 
within the errors with the values in Table |7] but the ratios ob- 
tained from the emission measure-independent method are more 
robust and have lower uncertainties. 

5. EQ Peg in the context of active l\1 dwarfs 

In order to put EQ Peg into the context of other similar late-type 
stars, we extracted five additional Chandra HETGS spectra of 
active M dwarfs for comparison with our results; our sample thus 
contains YY Gem, AU Mic, EV Lac, AD Leo, and Proxima Cen 
in addition to EQ Peg A and B. Most of these archival observa- 
tions have been included in many studies (e. g. lNess etal ■12003d. 
l2004t iTesta et akl 12007 ). the X-ray lightcurves and sp ectra of 
AU M ic, EV Lac, and Proxima Cen are shown by Test a et all 
(12004 . While the two stars in the binary EQ Peg are co-eval, we 
are now dealing with a range of different ages, rotation periods. 



and metallicity. In Table|9]we summarize for each star their spec- 
tral type, distance (from the SIMBAD database), age estimates 
and rotation periods from the literature together with observa- 
tion ID, exposure time, and X-ray luminosity obtained from in- 
tegration of the dispersed photons over the range 2-25 A. The 
bolo metric luminosity for Lv/Lhni has been calculated accord- 
ing to lKenvon & HartmannI ( 1995b . While the X-ray luminosity 
decreases along the sequence of spectral types, Lx/^^boi is in the 
saturation regime of w -3.3 for all stars except Proxima Cen. 

YY Gem is an MO. 5 companion of the Castor AB sys- 
tem. YY Gem itself is an eclipsing binary consisting of two al- 
most identical early M dwarfs with an orbital period of 0.814 
days at an inclination of 86.2, separated by 3.88 Rq. From 
eclipse mapping te c hniq ues based on an XMM-Newton obser- 
vation, iGUdeTeFal] (1200 1) found similar activity levels on both 
components, but concentrated in compact active regions at small 
scale heights, i. e. inter-binary loops as discussed for RS CVn 
system s are rather unlikely. This is confirmed by Stelzer et al.l 
(|2002|) . who modelled a flare in a simultaneous XMM-Newton 
and Chandra LETGS observation. Several strong flares occur 
also in the two HETGS observations discussed here. The young 
(w 10 Myrs) star AU Mic, spectral type Ml, has the highest 
X-ray luminosity within 10 pc to the Sun and has thus been ex- 
tensively observed in X-rays and in the EUV AU Mic shows 
a strong level of variability on all time sc ales and a high flare 
rate. From the HETGS dataset used here, iLinsky et al.l (l2002h 
computed an emission measure distribution with a peak tem- 
perature of log T X! 6.8 and subsolar abundances with a pro- 
nounced inverse FIP effect. Magee et al. (2003) found similar 
values for the abundances based on an XMM-Newton observa- 
tion. The Chandra HETGS observation of the M3.5 star EV Lac 
is part of a multiwavelengt h campaign and has been discussed 
in detail bv lOsten et al.l (2005, 2006). Several flares occurred to- 
wards the end of the observation, and lOsten et aP (12005 ) com- 
pared line ratios with different FIP dependencies during flares 
and quiescence and find a sligh t increase for low FlP/high FIP 
ratios during flares. lOsten et aTl (2006) derive a subsolar abun- 
dance level and an inverse FIP effect from an emission measure 
distribution of the quiescent st ate, with values consistent w ith an 
XMM observation analyzed by lRobrade & Schmit3(l2005l) . Both 
EV Lac and the M4.5 star AD Leo are young stars from the 
galactic di sk populatio n. AD Leo is another a well-known flare 
star (e. g. iHawley et al.l l2003l) . however, no larger flares have 
occurred during the Chandra HETGS observation discussed 
here. This dataset was analyzed by Maggio & Ness (2005), with 
the focus on density diagnostics. Maggio et al.i (i200 4) found a 
mild inverse FIP bias based on two Chandra LETGS observa- 
tions, which is roughly co nsistent with the results obtained by 
iRobrade & Schmita (l2005l) and from an XMM-Newton obser- 
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Table 9. Properties of active M dwarfs observed with Chandra HETGS. X-ray luminosities are computed from the total background- 
subtracted counts contained in the MEG spectra from 2-25 A. Values given in brackets correspond to the quiescent state. 



star 


spectral 


distance 


age 


PfOt 


ObslD 


exposure time 


logLx 


log 






type 


[pc] 


[ Myrs ] 


[d] 




[ks] 










YY Gem 


M0.5 


15.27 


370' 


0.814^ 


8504/8509 


136.0 


(81.5) 


29.63 


(29.43) 


-3.17 


(-3.37) 


AUMic 


Ml 


9.94 


122 


4.865^ 


17 


58.8 


(47.4) 


29.19 


(29.13) 


-3.48 


(-3.54) 


EQ Peg A 


M3.5 


6.25 


950^ 


(1.066^) 


8453/8484/8485/8486 


98.5 


(92.8) 


28.71 


(28.67) 


-3.19 


(-3.23) 


EVLac 


M3.5 


5.05 


453 


4.376* 


1885 


100.0 


(87.3) 


28.63 


(28.51) 


-3.10 


(-3.22) 


AD Leo 


M4.5 


4.69 


470^ 


27/149.10 


2570 


45.2 




28.54 




-3.36 




EQ Peg B 


M4.5 


6.25 


950^ 


(1.066') 


8453/8484/8485/8486 


98.5 




27.89 




-3.39 




Proxima Cen 


M5.5 


1.30 


5800-* 


83.5'° 


2388 


42.4 


(35.7) 


26.85 


(26.60) 


-3.96 


(-4.21) 



' [Torres & RibasI (|2002|) ^ IZuckerman et al.l (l200lh. ^ 
(l2007h . " iPettersen et alj ( ll992h . '^ ISpiesman & Hawlev 



Parsamvanl (fT99^, ^ lYildij iimlti . ^ iKronI <1952l) ^ [Torres et alj ( |1972|) , ' [Norton et al.] 
d 1 9861) . iMahmou j i 1 993h . ' ' [Benedict et al J d 1 998h 



vation. Ivan den Besse laar et al.l (l2003h analyzed the sani e data. 
While iRobra de & Sch mitti and Ivan den Besselaar et al.l deter- 
mined an overall subsolar abundance level, 'Maggio et al.' found 
values much higher than for the solar photosphere, however 
with a large uncertainty in the overall normalization. EQ Peg A, 
EV Lac, and AD Leo have comparable X-ray luminosities. 
On the basis of a common origin with the a Centauri system, 
Proxima Cen can be considered older than the Sun. With a spec- 
tral type of M5.5, it is also the latest star in our sample. Its 
X-ray luminosity is comparable to the Sun, more than two or- 
ders of magnitude lower than that of AU Mic. Nevertheless, 
Proxima Cen is also a well-known flare star. iGiidel et alJ 
(l2002ah and iGiidel et all (|2004 discussed a large-amplitude 
long-d urati on flare observed with XMM-Newton. Giidel et al. 
(12004 and iNordon & BehaJ (l2008h derived emission measure 
distributions for the diff'erent phases of the flare and found abun- 
dance ratios relative to iron to stay roughly constant, but with 
an overall tendency to an inverse FIP effect in the absolute level. 
Two flares are also included in the HETGS observation analyzed 
here. 

In order to obtain self-consistent results, we applied the same 
data reduction and analysis techniques as described above to 
the spectra of the other four sample stars. However, even with 
typical flare rates and amplitudes for active M dwarfs known, 
the occurrence of individual flares is randomly distributed, and 
each X-ray observation is a snapshot that is not necessarily rep- 
resentative of the typical flaring behavior of the observed star. 
We filtered the observation for quiescent-only emission and in- 
clude the quiescent luminosities in brackets in the last column of 
Tables g] 

We also computed the cumulative count spectra for the other 
five spectra as that shown in Fig.|3]for EQ Peg. In Fig.|9]we show 
these spectra for the quiescent states of the five stars as labeled in 
the legends, together with the well-exposed reference spectrum 
of Capella. For comparison, we included the cumulative spectra 
of EQ Peg A and B, bracketing the gray-shaded areas. In these 
plots, YY Gem, AU Mic and EV Lac are similar to the EQ Peg 
system, with YY Gem very similar to EQ Peg A, while AU Mic 
has stronger neon lines than EQ Peg, leading to larger steps in 
the cumulative distribution. AD Leo is more similar to EQ Peg 
B, and Proxima Cen is even softer. All seven M dwarfs dififer 
considerably from Capella. 

We applied our DEM-independent abundance determination 
using Eq. |6] with the coefficients listed in Table |6] to the total 
and quiescent-only spectra of AU Mic, EV Lac, AD Leo, and 
Proxima Cen and list the results in TablefTOl These numbers have 
to be compared to the relative abundances for EQ Peg A and B 
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Fig. 9. Cumulative distribution of counts in the quiescent spectra 
of the reference star Capella and other M dwarfs observed with 
the Chandra HETGS. The area between the curves for EQ Peg 
A and B is indicated in shades of gray for comparison. 

listed in Table [T] For comparison of our results with other abun- 
dance determinations for our sample stars, we tra n sformed the 
values determined by Robrade & Schmitt ("2005"), 'Giidel et aL 
(2001), Linskyetal. (2002), Magee et al. (2003), OstenetaL 
( 120061) . iMaggioet all (|2004), and iGudel etalj (|2()04|) to the ref- 
erence set of solar photospheric abundances froml AspIund et aH 
(l2005h and found good agreement within the errors. Only for 
neon we found a slight tendency of overestimating the abun- 
dance, as the blends to the neon lines may not always have been 
taken properly into account. 

For all seven stars, the Ne/O ratio is enhan ced by about a 
factor of 2, which is well-consistent with what iDrake & Testal 
(2005^) found for their sample of active stars. Si/Mg is also en- 
hanced by a factor of about 2, while Mg/Ne is depleted by factors 
from 3 (EQ Peg B and Proxima Cen) up to 10 (AU Mic). The 
values for O/N range from 0.3-1.0 times the solar level, partly 
with large uncertainties. The Ne/Fe ratio is clearly increased, by 
factors ranging from 4.5 (Proxima Cen) to 9 (AU Mic), while 
the Mg/Fe ratio is about solar. In Fig. [10] we give a graphical 
representation of the numbers listed in Tables|7]and[T0]from the 
quiescent data. 

O/N and Ne/O can be considered as high FlP/high FIP ratios, 
and within the errors, these ratios do not change with spectral 
type. No trends can also be identified for the low FlP/low FIP 
ratios Mg/Fe and Si/Mg. The low FlP/high FIP ratio Mg/Ne on 
the other hand clearly increases towards later spectral type, while 
the high FlP/low FIP ratio Ne/Fe decreases. 

Mg/Ne and Ne/Fe are sensitive to a possible FIP bias, and 
both are anomalously low and high, respectively, for all stars 
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Table 10. Abundance ratios of other M dwarfs observed with the Chandra HETGS. Values obtained for the quiescent state are given 
in brackets. 



ratio 


YYGem 


AUMic 




EV Lac 


AD Leo 


Proxima Cen 


O/N 


0.62±0.12 


(0.57+0.14) 


0.48+0.09 


(0.48= 


tO.ll) 


0.53+0.10 


(0.76+0.15) 


0.60= 


t0.13 


0.30+0.09 


(0.34+0.14) 


Ne/O 


1.97±0.12 


(2.10+0.18) 


2.36+0.19 


(2.31= 


t0.20) 


2.07+0.50 


(1.71+0.12) 


2.33= 


t0.19 


1.72+0.31 


(2.06+0.54) 


Mg/Ne 


0.18+0.02 


(0.16+0.02) 


0.11+0.01 


(0.12+0.02) 


0.24+0.02 


(0.22+0.02) 


0.21 = 


t0.02 


0.31+0.06 


(0.35+0.10) 


Si/Mg 


2.31+0.18 


(2.59+0.31) 


2.33+0.26 


(1.84= 


tO.26) 


2.47+0.22 


(2.25+0.25) 


2.46= 


t0.30 


2.48+0.28 


(2.50+1.28) 


S/Si 


1.27±0.18 


(1.46+0.33) 








1.82+0.33 












Ne/Fe 


7.08±0.33 


(7.34+0.46) 


8.85+0.57 


(9.03= 


tO.59) 


5.73+0.29 


(5.17+0.26) 


6.78= 


t0.41 


4.44+ 0.59 


(4.45+ 0.82) 


Mg/Fe 


1.23±0.19 


(1.24+0.24) 


0.84+0.18 


(1.03= 


tO.26) 


1.28+0.15 


(1.28+0.15) 


1.13= 


t0.19 


1.44+ 0.35 


(1.60+ 0.49) 



MO. 5 Ml M3.5 M3.5 M4.5 M4.5 M5.5 
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YY Gem AU Mic EQ Peg A EV Loc AD Leo EQ Peg B Prox Cen 

Fig. 10. Relative abundances as a function of spectral type de- 
rived from the seven M dwarfs as listed in TableslTjandfTOl 

in our sample, pointing towards an inverse FIP effect. This re- 
sult holds for both total and quiescent-only spectra. Since we 
only have relative abundances, it is unclear whether high-FIP 
elements are overabundant or low-FIP elements are underabun- 
dant. A low Mg/Ne ratio may be caused by an enhancement of 
neon, as supported by the increased Ne/O ratios, but also by a 
low Mg abundance, which would match the enhanced Si/Mg ra- 
tios, and so on. 

The trends observed for Ne/Fe and Mg/Ne are roughly inde- 
pendent of the occurrence of flares. This confirms previous find- 
ings that strength of the inverse FIP effect scales with the genera l 
activity level of a star (Gudel et al. 2002b; Telleschi et alfcoOSl), 
as activity (in terms of X-ray luminosity) decreases with increas- 
ing spectral type. While the behavior of the Ne/Fe and Mg/Ne 
ratios in our M dwarf sample implies an absolute decrease of 
neon with later spectral type, the constant Ne/O ratio does not 
fit that picture. This indicates that not only a single element like 
neon causes this effect, but that the inverse FIP effect in general 
diminishes for later M dwarfs. 

6. Discussion 

The behavior of M dwarf coronae shows striking changes with 
spectral type. The early M dwarfs up to a spectral type of 
^ M3-M4 observed so far are strong X-ray emitters, show- 
ing permanent flaring, hig h coronal temperatures a nd a pro- 
nounced inverse FIP effect dRobrade & Schmittl2005 >). X -ray lu- 
minosity, the level and peak temperature of the emission mea- 
sure, and the flaring rate decrease for later M dwarfs. While 
for earlier M subclasses, it is difficult - if not impossible - 



to distinguish between weak flaring and quiescence, M dwarfs 
of spectral type M5 or later show periods of quiescence inter- 
ruptedby discrete flares, and their coronal temperatures decrease 
dFuhrmeister et alJlIOOTi) . Very late M dwarfs (>M7) become 
more and more difficult to detect as quiescent X-ray sources, 
yet they can still produce enormous X-ray luminosities dur- 
ing flares. Abundance determinations are difficult for these ob- 
jects, but Giidel et al. (2004) and Fuhrmeister et al. (2007) found 
the inverse FIP effect to be less pronounced in the M5.5 stars 
Proxima Cen and CN Leo. With our systematic line ratio sur- 
vey, we can consolidate the trend of the inverse FIP effect in 
M dwarf coronae weakening with later spectral subclass. This is 
coupled to the decrease of the activity observed with later sub- 
class in general, while Lx/Lboi is similar and dominated by satu- 
ration effects for all our sample stars except for the latest object 
Proxima Cen. Apart from the precise strength of single emission 
lines that are mostly determined by individual abundances, high 
resolution X-ray spectra of M dwarfs have a very similar over- 
all shape. Only the continuum level decreases with later spectral 
type, as lower coronal temperatures reduce the amount of true 
continuum emission as well as the pseudo-continuum dominated 
by lines of highly ionized iron located at approximately 9-17 A. 
The reduction of the continuum level can easily be quantified 
from the cumulative spectra and serves as an indicator for the 
average coronal temperature. 

Both methods we used for our analysis, the cumulative spec- 
tra as well as the emission measure independent line ratios, can 
be applied to spectra with lower signal where only the strongest 
lines can be discerned. They are thus also suitable to determine 
coronal properties from the spectra of the fainter mid and late 
M dwarfs. Only very few grating spectra of such objects are 
available so far, and especially spectra of stars later than M6 
would help to confirm the trend that the difference of the Ne/Fe 
and Mg/Ne ratios to solar photospheric values decreases with 
later spectral type. However, since the Ne/O and Si/Mg ratios 
stay enhanced by about a factor of two, there is no indication 
of all abundance anomalies to diminish. So at least in compar- 
ison to the saturated coronae of active M dwarfs, it is the solar 
photosphere which appears anomalous; and the abundances of 
M dwarfs seem to "converge" with later spectral subclass, but to 
values differing from the solar photospheric level. 

The determination of reference sets of photospheric abun- 
dances of active M dwarfs is rather difficult, since reliable 
models of the stellar atmosphere including radiative transfer 
are required. The parameters involved in the modeling proce- 
dure, i. e. the stellar mass, effective temperature, log g, and for 
M dwarfs, the influence of the stellar chromosphere and pos- 
sibly the formation of dust, introduce large uncertainties in the 
derived abundances (Zboril & Byrne 1998). Additionally, only 
the abundances of heavier elements, often only that of iron, can 
be determined in the photosphere, while light elements like car- 
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bon, nitrogen and oxygen or noble gases like neon that dominate 
the coronal spectrum leave only weak or no signatures in the 
photospheric spectrum. 

It may be questionable to line up active M dwarfs exclusively 
based on spectral type, since other parameters, especially the 
age, certainly have a non-negligible influence, and age-activity 
relations are well-established over all wavelengths and in partic- 
ular f or coronal X-ray emission (Vaiana et al. 1992; Fa vata et al.l 
|1994|) . Our ultraluminous sample star AU Mic forms a proper 
motion pair with the M4/M4.5 binary AT Mic, both stars are 
members of the fi Pictoris rnoving group and only about 10-12 
Myrs old jZuckerman et all 1200 ll) . iRobrade & SchmittI (l2005h 
found AT Mic half an order of magnitude brighter in X-rays 
than our sample stars AD Leo, EV Lac and EQ Peg, which have 
similar spectral types. Stars with a common origin, like the two 
ZAMS stars AU Mic and AT Mic, or the co-eval components of 
the binary EQ Peg thus provide a testbed independent of such 
biases. The M3.5 and M4.5 components of EQ Peg differ by 
one spectral subclass, but the X-ray luminosity of the M4.5 sec- 
ondary is considerably lower, which we can attribute to a fairly 
reduced amount of coronal emission measure, but at tempera- 
tures similar to those found on the primary. Yet, in terms of 
Lx/Lboh both stars are in the saturation regime, and thus again 
comparable. Similar considerations apply to a larger sample of 
different ages, ranging from a few Myrs (AU Mic) to almost one 
Gyr (EQ Peg); the cumulative spectra of very young stars like 
AU Mic or EV Lac are very similar to the much older EQ Peg 
system. Proxima Cen is also by far the oldest object, and the 
one with the longest rotation period. It is the only star from our 
sample where Lx/L^oi deviates from the saturation level, yet it 
still fits the picture of the other saturated stars in terms of abun- 
dance trends. Here we probably see the well-established rela- 
tionships between X-ray luminosity and age or, related via spin- 
down with a ge, rotation period/rotational veloc ity and X-ray lu- 
minosity dPizzolato et al.ll2003t iReinersI [20071) or other activity 
indicators like Ha or Can emission at work. The age of the 
other sample stars is less than 1 Gyr and their rotation periods 
are typically only a few days. The sample of M dwarfs observed 
with high spectral resolution in X-rays is - by necessity - bi- 
ased towards the more luminous younger objects, and it would 
be crucial to examine how further old and inactive objects like 
e. g. Barnard's star fit into that picture. 

During quiescence, coronal densities on active stars typi- 
cally range between 10'"-10" cm for the low-temperature 
plasma (x 2 MK) traced by the O vii triplet jNess et al.ll2003ah . 
and between 10"-10'^ cm"^ for Neix (» 4 MK). Significantly 
higher values have s o far only been observed du ring huge flares 
(iGtidel et al.ll2002al: iLiefke" et aLllin preparation!) , but in general, 
coronal densities tend to increase during periods of increased 
activity and lower ampl i tude flaring (e. g.|M aggio & Ness 2005; 
iFuhrmeister et alj200'7h . lMitra-Kraev & Ness, (.2006) thus devel- 
oped a statistical approach that showed that the different ///- 
ratios observed during flares and quiescence on EV Lac are ac- 
tually highly significant. 

Flares can considerably change the X-ray properties of ac- 
tive M dwarf stars. The flare plasma often dominates the coronal 
emission, although the dimensions of the flaring loops seldomly 
reach the size of the star itself, and the flaring region typically 
covers only a small fraction of the stellar surface. The plasma 
temperature and emission measure can increase by more than 
an order of magnitude during the strongest flares, which sub- 
sequently results in an enhanced level of continuum emission. 
The bulk of line emission is dominated by high-temperature 
ions. Based on current flare models, the flare plasma is largely 



composed of evaporated chromospheric or photospheric mate- 
rial filling the loop. In this context, abundance gradients between 
the different layers of the atmosphere , like the FIP eff e ct ob- 
ser ved on the Sun, become important. iNordon & Behaii (l2007h 
an d lNordon&Behail (120081) investigated changes in the FIP bias 
during several larger flares on a sample of different active stars. 
For the very active stars that show an inverse FIP effect dur- 
ing quiescence, they found abundance ratios of the scheme high 
FlP/low FIP like Ne/Fe to decrease during flares, i. e. the inverse 
FIP effect turns into a solar-like pattern or even a FIP effect. This 
is in good agreement with the trends we find in our M dwarf 
sample. Thus abundance gradients, resulting in a different com- 
position of the evaporated flare plasma, seem to be a common 
feature, and in contradiction to the corona, the underlying abun- 
dances of the lower atmospheric layers in active M dwarfs ap- 
proach a more sol ar photospheric like pattern . This finding is 
also supported by Liefke et al.l din preparation!) , who found the 
coronal iron abundance of the M5.5 dwarf CN Leo to increase 
by a factor of two to the solar photospheric level during a large 
flare. 

To what extent flares do have to be considered as the "usual" 
behavior of a stellar corona? Early M dwarfs are in a state of per- 
manent flaring, and even during phases of apparent quiescence, 
the underlying basal coronal flux is difficult to disentangle from 
unresolved weaker flares. While larger flares can be separated, 
the characteristic flickering of early M dwarfs has to be taken 
as their typical "quasi-quiescent" behavior. The average coro- 
nal temperatures and densities of active stars will lower when 
the data is restricted to this quiescent state. However, we do not 
find that the flares included in the total datasets of our seven 
sample stars introduce large changes in the observed abundance 
ratios. Since the data do not include huge flares, this indicates 
a smooth transition between the quasi-quiescent state and in- 
termediate flares, which supports the assumption that the unre- 
solved flickering is indeed composed of smaller flares. When 
switching to intermediate and late M dwarfs and their single, 
discrete flares, it may on the other hand be possible to pin down 
the true basal coronal flux as observed on the Sun during qui- 
escence, which seems to approach solar coronal conditions for 
these less active stars. The smooth transition in X-ray luminos- 
ity, flare rate, coronal temperatures and abundances we observe 
in our sequence of very active early M dwarfs towards interme- 
diate and late M dwarfs, i. e., from stars with radiative interiors 
to fully convective interiors, suggests that the different dynamo 
mechanisms thought to operate in these stars to not lead to easily 
observable consequences in their coronal properties. 

7. Summary and conclusions 

We have investigated the coronal properties of the M3.5/M4.5 
EQ Peg binary system from their Chandra HETGS spectra. No 
large flares occur during our observations, and we find coronal 
densities ranging from 4 ■ 10'*^ cm""* to 3 ■ 10" cm"^ or consistent 
with the low-density limit from the O vii and Ne ix triplets, i. e. 
values typical for active stars. The ratio of X-ray luminosities 
is 6:1 for EQ Peg A and B, but in terms of Lx/ibob both stars 
are saturated. The differential emission measures of both compo- 
nents peak around 3 MK, and their abundances are similar and 
follow the inverse FIP effect. 

We compared all seven M dwarfs observed so far with the 
Chandra HETGS with two methods also suitable for spectra 
with low signal. The slope of the cumulative spectrum, which 
traces the continuum level and therefore the average coronal 
temperature, is very similar for spectral types MO. 5 to M4 and 
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then decreases. Emission measure-independent abundance ratios 
based on the line fluxes of strong lines with a similar temperature 
dependence confirm the existence of abundance anomalies com- 
pared to the solar photosphere for all M dwarfs in our sample. 
The ratios sensitive to a FIP bias, i. e. Mg/Ne and Ne/Fe, show 
a clear trend with increasing spectral type to approach the so- 
lar photospheric level, while ratios insensitive to a FIP bias like 
Si/Mg and Ne/O stay at a constant level. These trends seem to 
be independent of the age of the stars. 
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